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In this review we discuss the multifaceted problem of spin transport in hydrogenated graphene from a the¬ 
oretical perspective. The current experimental findings suggest that hydrogenation can either increase or de¬ 
crease spin lifetimes, which calls for clarification. We first discuss the spin-orbit coupling induced by local 
a — TT re-hybridization and sp® C-H defect formation together with the formation of a local magnetic moment. 
First-principles calculations of hydrogenated graphene unravel the strong interplay of spin-orbit and exchange 
couplings. The concept of magnetic scattering resonances, recently introduced ll|] is revisited by describing the 
local magnetism through the self-consistent Hubbard model in the mean field approximation in the dilute limit, 
while spin relaxation lengths and transport times are computed using an efficient real space order N wavepacket 
propagation method. Typical spin lifetimes on the order of 1 nanosecond are obtained for 1 ppm of hydrogen 
impurities (corresponding to transport time about 50 ps), and the scaling of spin lifetimes with impurity density 
is described by the Elliott-Yafet mechanism. This reinforces the statement that magnetism is the origin of the 
substantial spin polarization loss in the ultraclean graphene limit. 

PACS numbers: 72.80.Vp, 73.63.-b, 73.22.Pr, 72.15.Lh, 61.48.Gh 


Introduction.-Remaikahle electronic and transport prop¬ 
erties of graphene (see for instance 101) can be further tai¬ 
lored towards higher functionality through the use of chem¬ 
ical functionalization 10, irradiation (defect formation) 0 
or structural patterning (such as creating a nanomesh super¬ 
lattice) 0. Among the wealth of possible modifications, 
hydrogenation, ozonization or fluorination have demon¬ 
strated the tunability of graphene from a weakly disordered 
semimetal to a wide band-gap or Anderson insulator depend¬ 
ing on the nature and density of impurity atoms, varying 
typically from 0.001% to few percent iQ-lit]. Together with 
Raman spectroscopy which provides an estimate of defect 
density, analysis of the low-temperature conductance be¬ 
haviour and weak localization regime enable evaluation of 
the main transport length scales (such as the mean free path 
and localization length) ii. Surprisingly, strongly hydro¬ 
genated graphene, with anomalously large loffe-Regel ratio 
of 1/kpie ~ 240 still exhibits quantum Hall effect features 
in the high magnetic field regime d. 

Hydrogen defects are particularly interesting since, to¬ 
gether with a resonant scattering state created locally in 
space, the breaking of the sublattice symmetry entails the 
formation of a zero-energy mode and a local magnetic mo¬ 
ment on the order of 1 Bohr magneton for an isolated hydro¬ 
gen adatom d . Meanwhile, the coupling between different 


induced magnetic moments is either ferromagnetic or anti¬ 
ferromagnetic, depending on whether the H-adatoms corre¬ 
spond to the same or to different sublattices of the graphene 
lattice, respectively. Theoretical calculations have reported 
specific magnetoresistance si gna ls for specific long range 
magnetic ordering situations 11211 . which could be realized 
by substrate-induced chemical reactivity as proposed in Ref. 


by sul 

d. 


Graphene spintronics has attracted a lot of attention since 
the pioneering demonstration that spin could be efficiently 
injected and propagated over long distances at room tem¬ 
perature Q. This has opened an opportunity for the de¬ 
velopment of lateral spintronics lll5l - l2in . that would benefit 
from the unique features of graphene, such as a high mobil¬ 
ity and remarkable electronic and transport features like rel¬ 
ativistic energy dispersion, Klein tunneling phenomenon,... 
11231] . In Table 1, we compile the typical values obtained for 
charge and spin transport in a variety of graphene devices, 
from graphene supported on silicon oxide or boron nitride 
to suspended graphene. One observes that the variation of 
mobility varies by up to two orders of magnitude, whereas 
spin lifetime seems difficult to relate to the quality of the 
material. 


Understanding spin transport in hydrogenated graphene is 
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TABLE I: Charge/spin transport parameters in graphene on various substrates 


Substrate 

fj. {cm'‘/Vs) 

Dc (m‘‘s 

Ds (m^fi-®) 

Ts{ps) 

Xsfipm) 


2 X 10^ 

2 X 10"^ 

2 X 10-2 

~ 100 

1.5 - 2 

SiOzfH 

1 - 3 X 10® 

— 

1.3 X 10-2 

~ 500 

2.4 

Suspended|20ll 

3 X 10® 

0.1 -0.2 

0.05-0.1 

~ 150 

4.7 

SiCfn 

2 X 10® 

2 X 10-2 

4 X lO-"* 

2300 

0.5 - 1 

SiCfTil 

17 X 10® 

— 

— 

— 

> 100 

hBNfl^ 

4 X lO'^ 

0.05 

0.05 

200 

4.5 

hBN/G/hBN^ 

2.3 X 10^ 

0.03-0.1 

0.05 

3000 
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a challenging and important issue. The underlying goal is to 
use a low enough level of hydrogenation to preserve a siz¬ 
able transport signal while also inducing local or long range 
magnetic ordering (ideally a ferromagnetic ground state), 
and evaluate how spin diffusion is affected by interaction 
between itinerant spins and local magnetic moments. Ulti¬ 
mately, one could envision spin devices such as spin valves 
that do not use ferromagnetic materials to inject and detect 
spins, but rather utilize the varying signal response of mag¬ 
netized graphene to manipulate the spin degree of freedom 
and engineer logic functions, a long-standing quest of spin- 
tronics 01. 


From that perspective, the first experimental result show¬ 
ing some interaction between spin diffusion and magnetic 
moments produced by H adatoms was reported in 2012 by 
McCreary and coworkers ll25n . Low-temperature spin trans¬ 
port measurements on graphene spin valves (T=15K) were 
shown to exhibit a dip in the non-local spin signal as a 
function of the external magnetic field. This effect was 
tunable with hydrogen density, and was related to spin re¬ 
laxation spin relaxation by exchange coupling with para¬ 
magnetic moments. Spin lifetimes were estimated to be 
Tg ~ 400 — 500ps in absence of hydrogen, and were slightly 
enhanced upon increasing the H-adatom density up to nn ~ 
0.1%), and accompanied by a reduction of charge mobility 
by one order of magnitude to /i ^ 500cm^/Vs. 


A subsequent experiment also reported a surprising in¬ 
crease of spin lifetime with hydrogen density, with values 
enhanced by a factor of two after hydrogenation of the sam¬ 
ples for uh ~ 0.02% ll26ll . The upscaling of Tg with hy¬ 
drogen density suggests a Dyakonov-Perel relaxation mech¬ 
anism ll27n . In these experiments, the effect of magnetic 
moments is proposed to be strong enough to counteract the 
expected increase of spin-orbit coupling, which would be 
expected to produce shorter spin lifetimes with increasing 
H-adatom density. 


In 2013, J. Balakrishnan and coworkers lE^ reported ex¬ 
perimental evidence for a room temperature spin Hall ef¬ 
fect (SHE) in weakly hydrogenated graphene, with non¬ 
local spin signals up to lOOfl (orders of magnitude larger 
than in metals). The non-local SHE revealed by Larmor 
spin-precession measurements was assigned to a colossal 
enhancement of the spin-orbit interaction induced by H- 


adatoms for density in the range uh ~ 0.01% — 0.05%, 
mobilities /r ~ 14.000 — 900cm^/Vs) and spin lifetime 
on the order of Tg ^ lOOps. The spin-orbit interaction 
was estimated to be about 2.5 meV, one order of magni¬ 
tude higher than ab initio calculations 1^ . The observa¬ 
tion of the large SHE-signal was exclusively assigned to an 
enhancement of spin-orbit coupling by H-adatoms, which 
limits spin lifetimes. It is however not clear whether the pre¬ 
viously observed contribution of magnetic moments in spin 
transport |25, 2^ remains marginal or not for the forma¬ 
tion and strength of the SHE signal, and which mechanism 
(Dyakonov-Perel lE^ or Elliott-Yafet |0]) can explain the 
variety of conflicting experimental data. 

In this context, it is of prime importance to clarify the spe¬ 
cific impact of H-adatoms on spin relaxation in graphene, 
and to quantify the relative contribution of spin-orbit cou¬ 
pling and magnetic effects. A first fundamental advance in 
this direction has been made by Kochan, Gmitra and Fabian 
ifUl who introduced a new relaxation mechanism in hydro¬ 
genated graphene driven by resonant scattering by magnetic 
impurities. By neglecting the spin-orbit coupling effect, spin 
lifetimes of Tg ^ 150ps were estimated for 1 ppm Hydrogen 
at room temperature and for large contribution of electron- 
hole puddles. 


LOCAL MAGNETIC MOMENTS AND SPIN-ORBIT 
COUPLING PROM PIRST-PRINCIPLES 


Here we remind the main ii^redients and predictions of 
the model introduced in Ref.lUj], and present new results 
obtained with a more generalized microscopic approach to 
study spin dynamics and relaxation effects induced by mag¬ 
netic moments. In the latter model, we describe the local 
moment formation using the self-consistent Hubbard model 
on hydrogenated graphene lll2ll . This model is shown to 
correctly reproduce ab initio results for spin splittings ob¬ 
tained for small supercells, whereas it provides, by a scaling 
analysis, a more extended spreading of local magnetic mo¬ 
ments around the H ii™urity when compared with the an¬ 
alytic model used in yj]. Next, the Hubbard tight-binding 
parameters are implemented into a real space wavepacket 
propagation method which gives direct access to spin dy- 
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namics and spin relaxation effects 0. 


First-principles calculations of local magnetism and spin-orbit 
coupling effects in hydrogenated graphene 

When a hydrogen atom is chemisorbed, graphene locally 
undergoes a structural deformation, breaking the sp^ sym¬ 
metry in favor of an sp^-like hybridization. After ab ini¬ 
tio structural optimization, the hydrogenated carbon atom is 
found to be slightly displaced out of the plane 0.5A), 
forming a C-H bond of length 1.13A. As measured experi¬ 
mentally, the (T — TT rehybridization induced by the hydrogen 
adatom comes with the formation of a local magnetic mo¬ 
ment of the order of 1 pB ■ Recently, spin-unpolarized first- 
principles calculations suggested that the local sp^ distortion 
and its resulting pseudospin inversion asymmetry induce a 
giant enhancement of the spin-orbit coupling (SOC) 0 . In 
this section, spin-polarized first-principles calculations are 
used to investigate the interplay between the local magnetic 
moment and the enhanced spin-orbit interaction induced by 
hydrogen adatoms. Using unconstrained spinors to represent 
the one-particle wave-functions, the respective contributions 
to the energy bands splitting of both SOC and electronic ex¬ 
change are computed. The spin textures associated with the 
low energy electronic spectrum are also investigated. 

Hydrogenated graphene in the dilute limit is represented 
by a large supercell (5x5x1) containing a single hydrogen 
defect, leading to a hydrogen concentration of ~ 2%. The 
out-of-plane dimension of the cell ensures a distance > ISA 
between neighboring graphene planes in order to avoid inter¬ 
action between periodically repeated images. The one parti¬ 
cle Hamiltonian is computed within the framework of non- 
collinear spin-polarized density functional theory as imple¬ 
mented in the Vienna Ab initio Simulation Package (VASP) 
0. The projector augmented wave method is used to ex¬ 
pand the one-particle wave-functions up to an energy cut¬ 
off of 600eV i^ . The eigenstates of the self-consistent 
Hamiltonian are populated according to an electronic tem¬ 
perature of 300K. Electronic exchange and correlation are 
treated within the generalized gradient approximation by 
means of the PBE functional 13311 . Integrals over the Bril- 
louin zone are performed using a 6x6x1 Monkhorst-Pack 
grid of k-points. The geometry is fully optimized until re¬ 
maining atomic forces and stresses are lower than O.OleV/A 
and O.OleV/A^ respectively. Various magnetization axes, ei¬ 
ther parallel or perpendicular to the graphene plane, have 
been considered. The different configurations are energeti¬ 
cally equivalent, i.e. the computed total energies vary by less 
than 1 peV. In what follows, we focus on the most symmet¬ 
rical and experimentally relevant configuration where the 
magnetization axis is chosen as the direction perpendicular 
to the graphene plane. 

The computed low-energy band structure is depicted in 



(a) 



r M K r 

(b) 
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FIG. 1: (a,b) Ab initio band structures of hydrogenated graphene 
in the dilute limit 2%) including spin-orbit coupling calcu¬ 
lated within both (a) spin-polarized and (b) spin-unpolarized frame¬ 
works. The color scale is associated with the in-plane component 
of the spin expectation values. Hence, the blue and red extrema of 
the color scale correspond to eigenstates whose spin is parallel and 
perpendicular to the graphene plane. Black dots are determined by 
symmetry. (c,d) In plane projection of the spin texture of the high¬ 
est valence band in both the spin-polarized (c) and -unpolarized (d) 
frameworks. The arrows indicate the in-plane projections of the 
spin directions (rescaled for eye convenience). The color scale in¬ 
dicates the actual values of the in-plane projections. Red and blue 
arrows respectively correspond to spin momenta that are parallel 
and perpendicular to the graphene plane. 


EigH] (a). The hydrogen impurity opens an energy gap 
of 1.24eV and produces two nearly-flat bands around the 
Eermi energy, accounting for the computed total magnetic 
moment of 0.99/iB. In this configuration, the spin degen¬ 
eracy is broken and the energy splitting of the bands origi¬ 
nates from both electronic exchange and spin-orbit interac¬ 
tion. Eor the sake of comparison, the low energy spectrum 
computed within a spin-unpolarized framework is also re¬ 
ported in Eig[T] (b). We note that spin-unpolarized means 
that the density matrix is constrained to be diagonal with 
same values for both spin-directions. In this case, the energy 
splitting of the bands only originates from the spin orbit in¬ 
teraction. At the scale of EigUlb), the splitting is not visible 
to the eye. 

The computed energy splittings corresponding to the low- 
energy bands are reported in Eigl2](a,c). The energy splitting 
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computed within the spin unpolarized framework [see Fig|2] 
(c)] is in perfect agreement with previous estimates lE^ . In 
order to clarify the impact of the spin-orbit interaction, the 
SOC-induced splitting here defined as the difference 

between the energy splittings computed with and without the 
spin-orbit interaction is also reported in Fig|2](b). is 

shown to be at least two orders of magnitude smaller than the 
energy splitting associated with the stabilization of the lo¬ 
cal magnetic moment by the electronic exchange (Figl2](a)). 
This highlights the prevailing role of electronic exchange in 
ruling the low-energy spectrum in the presence of local mag¬ 
netic moments. The comparison of Figl2](b) and (c) further 
emphasizes the interplay between spin-orbit and exchange 
couplings. 

The reduction of SOC-induced spin splitting due to ex¬ 
change coupling is explained by considering a simple ef¬ 
fective Hamiltonian. Consider electrons with momentum k 
along X and energy Sk, experiencing Rashba SOC a which 
is in general momentum dependent. If the electrons also feel 
an exchange coupling A with a magnetic moment oriented 
along z, as is our case, the Hamiltonian is 

H = Ek + Affz + a{k)ay. (1) 

The energy splitting of the eigenvalues of this operator is 

2\/A2 -Pa2. (2) 

In the absence of exchange, the splitting is linear in a, since 

Rashba SOC splits two degenerate levels. In the presence 

of exchange, and for realistic A ^ a, as in our case, the 

Rashba contribution to the splitting is only second order in 
2 

a (^)- This rough estimate gives for SOC of 100 peV and 
A of 0.1 eV, a Rashba coupling contribution to the splitting 
of approximately 0.1 peV, significantly below what is ob¬ 
served in Fig 12] (b). It is then likely that the difference in the 
calculated exchange coupling with and without SOC is due 
to the Rashba dependence of the exchange itself. This could 
be linear with a, with a numerical factor. The understand¬ 
ing of this dependence in terms of simple phenomenological 
models is still elusive and is currently the object of further 
investigation. 

The spin texture of the valence band is illustrated in FiglU 
(c). When the electronic density matrix is constrained to be 
spin-unpolarized, the frustration of the electronic exchange 
’’artificially” leads to a highly non-collinear spin texture en¬ 
tirely ruled by SOC (see Fig[T](d)). On the contrary, in the 
unconstrained framework, the stabilization of the magnetic 
moment dominates the energetics of the low energy spec¬ 
trum and the spin texture is found to be nearly collinear with 
the spin-orbit interaction only accounting for ^ 0.1% of the 
deviation from collinearity (FiglT](c)). While the local rehy¬ 
bridization of graphene induces an enhancement of the spin- 
orbit interaction, the magnetic properties of graphene as well 
as its spin texture are quasi-exclusively ruled by the elec¬ 
tronic exchange in the presence of hydrogen-induced local 
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FIG. 2: Energy splittings of the conduction (red curves), impu¬ 
rity (blue curves) and valence (green curves) bands along high- 
symmetry lines in reciprocal space, (a) Total energy splitting within 
the spin -polarized framework, (b) Contribution of SOC to the en¬ 
ergy splitting within the spin polarized framework (i.e. difference 
between the energy splittings computed respectively with and with¬ 
out accounting for the contribution of the spin-orbit interaction to 
Hamiltonian), (c) Total energy splitting within the spin-unpolarized 
framework (i.e. only arising from SOC as the electronic density is 
constrained to be spin-up polarized). 


magnetic moments. Hence, the contribution of the spin-orbit 
interaction to the Hamiltonian is neglected in the following 
discussion on the spin lifetimes in hydrogenated graphene. 


ONE-ORBITAL MEAN-FIELD HUBBARD 
APPROXIMATION 


The electronic structure of graphene is modeled by a nearest- 
neighbor tight-binding model with a single p^-orbital per 
site. When a H atom is adsorbed on top of a carbon (C) 
atom, the sp^-symmetry is locally broken, and the electron 
from the C Pz orbital is removed from the tt bands to form 
a (T bond with the H atom (Figl3](a)). Figure Ob) shows a 
5x5 supercell where the absence of a p^-orbital in the cen¬ 
ter of the green region (yellow site) represents hydrogen ad¬ 
sorption. To remove the Pz orbital, we use a sufficiently 
large on-site potential V°° ~ 10"^ eV. To properly describe 
magnetism in hydrogenated graphene, we introduce on-site 
Coulomb repulsion between electrons with opposite spins by 
means of the Hubbard model in its mean-field approximation 

= -70 + UY^ , 

i 

(3) 

where t is the first-neighbors hopping term, cl^ (cj,cr) is the 
creation (annihilation) operator in the lattice site i (j) with 
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spin a, U is the on-site Coulomb repulsion, and 
are the converged expectation values of the occupation num¬ 
bers for spin-down and spin-up electrons, respectively. The 
ratio t// 7 o = 1 has been chosen to accurately reproduce 
first-principles calculations. 




FIG. 4; Evolution of magnetic moments on first neighbor atom for 
increasing distance between impurities {(Ih-h)- The value con¬ 
verges apparently for (Ih-h > 33.74 A where long range mag¬ 
netic interactions start to vanish precluding the formation of FM or 
AFM ordering and leading to isolated paramagnetic defects. On 
the right side of the figure, we show the magnetic moments m; 
corresponding to each lattice site for different supercell sizes. 


FIG. 3: (a) Schematic representation of the sp^ hybridization 
breakdown upon hydrogen adsorption, equivalent to the creation 
of a single vacancy, (b) Single-vacancy 5x5 supercell used in 
the mean-field Hubbard calculations with periodic boundary con¬ 
ditions. The pz orbital at the center of the vacancy (yellow site 
inside the green region) is removed by adding a large on-site po¬ 
tential (17°“) to simulate hydrogen adsorption. 


Hydrogen defects in graphene may induce sublattice sym¬ 
metry breaking leading to the appearance of zero-energy 
modes in the density of states (DOS), as predicted by Inui 
et al. |4ffl, which are mainly localized around the impu¬ 
rities II 35 I [ 3 ^. These zero-energy states become spin- 
polarized upon switching on the Coulomb repulsion, lead¬ 
ing to semi-local S = 1/2 magnetic moments with a 
stag gered sp in density primarily localized on one sublat¬ 
tice il ll I 371 13811 . For a finite concentration of defects, the 


long range ordering of magnetic moments is dictated by 
the type of sublattice functionalization, being co-polarized 
or ferromagnetic (FM) for the same sublattice and counter- 
polarized or antiferromagnetic (AFM) otherwise. The total 
spin S of the macroscopic ground state is given by the ex¬ 
cess of magnetic moments on one chosen sublattice 1^ . 
although S' = 0 is the most likely value on simple statistical 
grounds (equal H occupation of both sublattices). Here we 
consider the dilute limit so that long range magnetic states 
are neglected and we get local paramagnetic (PM) impuri¬ 
ties. 

In order to reach a pure paramagnetic state, we have stud¬ 
ied the evolution of the local magnetic moments rm = 
1/2 — rii^i) with the varying supercell size (FiglU). In 

particular, Figl4](left panel) shows the values of the magnetic 
moment only at nearest neighbor sites of the defect m„„. 


As the supercell size is increased and the defects move away 
from each other, decreases quickly and converge for a 
supercell size 14x14, which corresponds to 0.25 % of ad¬ 
sorbed H and dn-H = 33.74 A. This decay is related to the 
overlap between wavefunctions corresponding to neighbor¬ 
ing defects. For strongly overlapping states (large concen¬ 
tration of H), long-range interactions induce a FM state be¬ 
tween magnetic moments located in neighboring cells which 
increases the values of rm. When the overlap is very small 
(diluted limit), long-range interactions have no influence on 
the local magnetism leading to almost constant values 
around the defect (PM state). 

Figure |5ja) shows the band structure of a single hydro¬ 
genated 28x28 graphene supercell, corresponding to 0.06% 
of adsorbed H, calculated using the mean-field Hubbard 
model (Eq. [^l. The splitting (Ag) of the mid-gap state 
formed during H adsorption is plotted in Fig. Qb) for dif¬ 
ferent supercell sizes at different fc-points. At the K point 
Ag decays as oc while at the F and M-points Ag oc 

These results confirm the existence of very small 
splittings at tiny concentrations of H (of the order of Ippm). 
In order to compare the values of the energy splitting of 
spin polarized bands obtained using the mean-held Hubbard 
model with published DFT results ll]], we have plotted Ag 
along M-K-F path in the inset of Fig. |3b). Bottom, mid¬ 
dle and top panels correspond to the valence, mid-gap and 
conduction bands respectively. Although DFT results for a 
5x5 supercell (dashed lines) show a slightly higher splitting 
of the spin polarized bands with respect to our results, this 
methodology nicely capture the main physics of magnetic 
resonances induced by H adatoms in graphene. 
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dH-H(A) 


FIG. 5: (a) Band structure corresponding to 0.06% of H adsorbed 
(28x28 supercell). The spin degeneracy is broken at the Fermi en¬ 
ergy (E = 0) due to the formation of local magnetic fields around 
impurities, (b) Evolution of As at three symmetry points K, F 
and M, for increasing values of the distance between impurities 
(d/r-ff). The inset shows the evolution of As along the M-K-F 
path in the hexagonal Brillouin zone of the valence (bottom), impu¬ 
rity (center) and conduction (top) bands. Dashed lines correspond 
to DFT results. 


SPIN DYNAMICS AND RELAXATION 


from the ab-initio calculation in Refs. Qlil —the hydro¬ 
gen on-site energy eh = 0.16 eV and graphene-hydrogen 
hybridization T = 7.5 eV. The exchange term —Jsh-S de¬ 
scribes an effective interaction of the local magnetic moment 
S induced on the hydrogen site with an itinerant electron 
spin Sh when this hops into the hydrogen \h) = /i+|0)-level 
that hybridizes with the graphene host. Both, S and Sh are 
vectors of spin i-operators (in our definition without con¬ 
ventional I factor) and J is a constant with the dimension 
of energy. For the exchange we take J = — 0.4eV since 
this value is consistent with a more detailed parametrization 
of the magnetic impurity model as discussed in Ref. HI]. In 
fact the precise value and the sign of J are not important as 
long as I J| is greater that the orbital resonance width (for the 
case of hydrogen the resonance width is 4 meV). In the inde¬ 
pendent electron-impurity picture (we do not discuss Kondo 
physics), the exchange term — Js/j • S can be diagonalized 
introducing the singlet (i = 0) and triplet (£ = 1) composite 
spin states |£, me) {me = —£, ...,£). Transforming the ini¬ 
tial Hamiltonian TLs, Eq. (jUi, to the new spin basis we arrive 
at 

Us = 0 ^ [e/i - J[‘^£ - 3)] \£,mt){£,me\ 

(5) 

+ T 0 E \£,me){£,me \. 


We now investigate the spin dynamics and relaxation phe¬ 
nomena in hydrogenated graphene by comparing two com¬ 
plementary theoretical approaches and contrasting our re¬ 
sults with state-of-the-art experimental data. 


Single impurity limit 


First, we consider the magnetic scattering problem in the 
single impurity limit, which can be simplified by consider¬ 
ing only the spin-flip processes stemming from the exchange 
coupling on the resonant scatterer (H atom) site. This was 
done in details in Ref. Ill] on the basis of the Hamiltonian 
TL — 77o +T~Ls which involves the usual -orbital Hamil¬ 
tonian for the pristine graphene T-Iq = — 7 o S<ij) 
together with the term T-Ls describing graphene-hydrogen 
chemisorption including the interaction of electron spin Sh 
and impurity moment S (exchange coupling) 


(7 

-Jsh-S 


T 


ht 


<7 H 


+ c: 


Ch,ct' 


(4) 


Here ) and c+ ^ (c^ ^) are fermionic operators creat¬ 

ing (annihilating) electron with spin cr on the hydrogen and 
graphene carbon site m, respectively. The orbital chemisorp¬ 
tion parameters entering the Hamiltonian TLs are extracted 


Down-folding the hydrogen |/i)-state by means of the 
Lowdin transformation we get for each spin component 
£, me an independent effective delta-function problem 

TifiE) = 0 \£,me){£,me\, 


located on the hydrogenated carbon site Ch with the energy 
dependent coupling 


Vr(E) 


T2 

E-eh + J{U - 3) ■ 


(7) 


From here it is straightforward to compute T-matrix ele¬ 
ments for Bloch states |k) and ]/«') of the unperturbed 
graphene (to shorthand the notation n comprises band in¬ 
dex and crystal momentum). Assuming the Bloch states are 
normalized to graphene unit cell the result is as follows 


,m^r{^E) ^ .y^ (^E)G (F) ^ 

where Ye{E) is given by Eq. (j?]) and Go{E) is Green’s func¬ 
tion per site and spin for the unperturbed graphene, i.e., 




E^ 




-,^]^q(W-\E\). (9) 


The above Green’s function is valid near the graphene charge 
neutrality point in the energy window from -1 eV to 1 eV, 
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where the linearized bandwidth W = VV^7o — 6 eV. 
For practical reasons we need relaxation rates 1 {E) 

that characterize spin-conserving (cr = a') and spin-flipping 
(cr ^ a') processes in graphene at the given Fermi energy E 
in presence of magnetic active impurities. For that we take 
into account all scattering processes |k, cr, S) —>■ |k', cr', E') 
at the given energy E = E{k) with the requested incom¬ 
ing and outgoing electron (hole) spins cr and cr' and allowed 
impurity spins E and E' (a charge carrier in graphene flips 
its spin only if the magnetic moment does the same to con¬ 
serve the total angular momentum). To get 1/tuu'{E) we 
start with the transition rate trace out the 

E-spin degrees of freedom corresponding to magnetic mo¬ 
ment S induced on the hydrogen site, then we integrate over 
all outgoing momenta |k') and finally we average the result 
over all incoming states |/«) at the given Fermi energy E. 
This can be done because the angular dependence of the T- 
matrix is trivial, it does not depend on the angle of k neither 
of k'. Assuming the distribution of magnetic moments is 
dilute and they are in average spin-unpolarized (e.g. by in¬ 
teraction with phonons) the result is as follows 


perature T then becomes 

, E[-5..F(efe,F;,T)]l/r.(efe) 

-= — - ( 12 ) 

k 

where F{ek, E,T) is the equilibrium Fermi-Dirac distri¬ 
bution at temperature T and Fermi energy E. Experi¬ 
ments show that charged impurities and electron-hole pud¬ 
dles within the sample still further affect system’s Fermi en¬ 
ergy. To account for all that effects one usually convolves 
1/ts{E, T) by the Gaussian kernel with standard deviation 
(Tb, 

oo 

1 ^ 1 f , 1 f 

Ts{E,T) 2^/TTcrb J ^Ts{e,T)^^^\ 2al )' 

— OO 

(13) 

Depending on the sample quality the typical value of crt, to 
be taken for best fits is about 70-110meV (see discussion 
next section). 


^ ^27r,. cpw r Vi(g) _ Vo(g) 

Taa'{E) ~ ^ ^ [ i_Vi(£;)Go(E) ’ l-Vo{E)Go{E) _ ' 

( 10 ) 

Here p is the concentration of hydrogen impurities per car¬ 
bon atom, Vf){E) = — F ImGo(£') is the graphene density 
of states per atom and spin, Vg{E) for f = 0,1 is given by 
Eq. 0 and function f^a' [a^, y] is defined by 

1 I I 2 1 I I 2 

f<ya'[x,y\ =-5aa'\x\ + -|a; + (ct • cr') y| . (11) 

The function f^a' originates from the decomposition of the 
actual spin states |tT, E) and |(t',E') with respect to the 
singlet-triplet basis |f, m^) after the tracing out the spin 
states of the induced magnetic moment. The symbol a ■ a' 
entering its definition equals 1 (-1) for the parallel (antipar¬ 
allel) spin alignments and | | stands for the absolute value. 

Knowing the partial rates XjTfjfjiiE) we can define the 
spin-relaxation rate \/Ta{E) = 1/t\,^{E) -\- \/t^^[E) and 
momentum-relaxation rate 1 /t(E) = l/T-|-|-(iJ)-|-l/T'|'j^(£’) 
at zero temperature as functions of Eermi energy. Erom 
Eq. (ITOl i we immediately see that both l/r^ and l/r should 
go to zero at the charge neutrality point since there are no 
states that can participate in scattering, (0) = 0- Secondly, 
the denominators 1 — Yo{E)Go{E) and 1 — Yi{E)Gq{E) 
that enter the function /o-o-' can minimize at certain energies 
Eq and Ei and this would manifests as two sharp peaks (sin¬ 
glet and triplet one) in the spin and momentum relaxation 
rates. Both mentioned features are clearly seen in Eig.|7] 

To account for finite temperature effects one should ther¬ 
mally broaden I/ts implementing rate equations for the 
graphene charge-carriers that obey Eermi-Dirac statistics. 
The concise formula for the spin-relaxation rate at finite tem- 


Tight-binding Model of hydrogenated graphene 

In our second approach to spin relaxation, we study the 
spin dynamics using the following tight-binding Hamilto¬ 
nian, 


'HSr = ^h'^Knhm+T^h'i;^Ci (14) 

{mi) 

+ JjcfazCi 

i 

which includes the self-consistent Hubbard terms describ¬ 
ing isolated magnetic moments (dilute limit), and a ran¬ 
dom distribution of H-adatoms in the ppm range. The long 
range nature of magnetic moment induced by hydrogen is 
involved in this Hamiltonian by considering up to ninth- 
nearest-neighbor exchange coupling term J^. Spin-orbit 
coupling is found to yield negligible corrections to the re¬ 
sults, so it is neglected in the following. The spin dynamics 
are investigated by computing the time-dependence of the 
spin polarization defined by 1301] 


S{E,t) 


Tr 

S{E-H)s{t) 

Tr 

5{E - H) 



where s(t) = e~^se~^. Using random phase states to 
perform the trace efficiently, we get 

S(i5, t) = - H)s + sSiE - H)m)) 




















and the initial wavepacket can be prepared in an arbitrary 
spin polarization as 



In what follows, the wavepackets are prepared (at t = 0) 
with an in-plane spin polarization and their time-dependent 
propagation and Sx{E, t) are evaluated applying the evolu¬ 
tion operator (Schrodinger equation) to the wavepackets . In 
absence of spin-orbit interaction, the spin dynamics is influ¬ 
enced by the existence of local magnetism, with H-related 
magnetic moments pointing out-of-plane. 

Fig|6] shows the time-evolution of the spin polarization at 
three different energies, namely the Dirac point, an energy 
close to the expected resonance and some high energy value 
(see inset). The curves first exhibit an sudden drop (espe¬ 
cially for the Dirac point) followed by an exponential decay, 
which dictates the values of the spin relaxation times Ts us¬ 
ing Sx{t) = (fitting the numerical results from 

fo = 75ps). The initial fast decay of Sxit) is understood 
as follows. We study the propagation of wavepackets which 
(for computational efficiency) are at time f = 0 in a random 
phase state, that is extended through the whole system. For 
smooth disorder, this method well captures the main trans¬ 
port length scales as described in Ref.||2l. Here however, 
the hydrogen adatoms (or vacancies) produce strong mag¬ 
netic moments localized around the impurity. This intro¬ 
duces some transient decay of the initial spin polarization 
which is not representative of the long time exponential be¬ 
havior. 


Spin relaxation times in hydrogenated graphene 

In Figl?] the inverse spin relaxation times obtained for 
the two models are superimposed. Within the single im¬ 
purity approximation, two sharp magnetic resonances are 
clearly formed (for singlet and triplet states), at which 1 /ts 
is maximum at zero temperature (green dashed lines). This 
approximation of Tg however diverges close to the Dirac 
point, which disagrees with most low-temperature exper¬ 
imental data. Room temperature spin-relaxation experi¬ 
mental data can however be reproduced using thermal and 
charge-puddles broadenings. A choice of 0.36 ppm of H- 
impurities together with a finite temperature broadening 
(300 K) and charge density fluctuations of 110 meV are 
necessary (blue dashed line) to reproduce experimental data 
(black filled circles) obtained for high-quality graphene spin- 
valve devices. lE^ . 

The inverse spin lifetimes obtained by the real space spin 
propagation method are shown in Fig|2]for a H-density of 1 
ppm (black solid line), 2 ppm (red solid line) and 15 ppm 
(blue solid line). The presence of the two resonances is con¬ 
firmed by this calculation. At the resonances, the inverse 



FIG. 6: (Main frame) time-dependent polarization for initial in¬ 
plane spin polarized wavepackets at selected energies (see inset). 
Numerical fits using an exponential damping are shown (dashed 
lines) on top of numerical simulations start from elapsed propaga¬ 
tion time t > 75ps. Inset: total density of state for the case of Ippm 
of hydrogen on graphene. 


spin lifetime is greater than in the single impurity approx¬ 
imation, which is likely due to the different natures of the 
adatom models. In the single-impurity adatom model the the 
magnetic moment is strongly localized and the resonance is 
sharp, while in the Hubbard model we have an effective de¬ 
fect model with more delocalized magnetism is real space. 
This broadening reduces the “residency time ” of the elec¬ 
trons on the magnetic virtual bound state (as pictured in the 
inset of Fig|2l-inset), and thus yields longer spin lifetimes. 
Indeed for 0.36 ppm of H-impurities, the two models differ 
to three orders of magnitude in spin lifetime. By further con¬ 
trasting the real-space spin results to the experimental data 
(black dots in Fig|7]i, we extrapolate that a density of 10-20 
ppm of H-adatoms can reproduce the range of experimental 
values of the graphene samples measured in Ref. 12611 . 

To discuss the relaxation mechanism, we compute the mo¬ 
mentum relaxation time Tp first from the sum of spin con¬ 
serving rate and spin-flip rate, which both can be calcu¬ 
lated by transforming the singlet and triplet T-matrix am¬ 
plitudes via composite spin states of electron and impu¬ 
rity in. The second approach uses the full Hamiltonian 
ids and the time dependence of the diffusion coefficient 
D{E,t) H. For the hydrogen coverage as low as 0.02%, 
we still can observe the saturation of the diffusion coefficient 
D{E, t) —> Djnax{E) , allowing us to extract Tp as 


Tp{E) 


Etnax (^E'j 
4u2(f;) 


(15) 


where v{E) is the pristine graphene velocity. We use the 
Fermi golden rule and the expected scaling of Tp{E) ~ 
1 / uh, to extrapolate the values for each ppm concentration 









E (meV) 


FIG. 7: (Main frame) Energy-dependence of Ts derived from the 
T-matrix of the single impurity model at zero temperature (green 
dashed lines); and at 300 K broadened hy puddles with energy fluc¬ 
tuations of 110 meV (blue dashed lines). Same quantity obtained 
from the real space spin propagation method for 1 ppm (solid black 
line), 2 ppm (solid red line) and 15 ppm (solid blue line) of H 
adatoms. Inset: schematic of the magnetic resonance process at 
the origin of enhanced spin relaxation. 


(see Figl8]l. The result is given in FigH] (green solid lines), 
and shows a weak energy dependence, with slight increase 
close to the Dirac point, with values in the range of 60-70 ps 
for 1 ppm of hydrogen impurities. In FiglS^inset), for the 
two models are reported showing some discrepancy close to 
the Dirac point, as expected from the approximations made. 
The values match well at high enough energy. 

We observe that Tg ^ lOr^, while the scaling of with 
impurity density clearly manifests an Elliott-Yafet spin re¬ 
laxation mechanism, as predicted for such type of impuri¬ 
ties 114ill . In some experiments, Ts is found to increase with 
hydrogen concentration, suggesting differently a Dyakonov- 
Perel relaxation mechanism lE^ . whereas others report on 
the opposite trend associated to the Elliott-Yafet mechanism 
0281] . The origin of such inconsistency remains obscure at 
that point but could stem (in particular) from the segrega¬ 
tion of H-adatoms and the alteration of fundamental effects 
induced by magnetic resonances. 

In conclusion, we investigated the impact of hydro¬ 
gen adatoms on charge and spin transport in transport in 
graphene in the dilute limit (down to the ppm limit), for 
which no long range magnetic ordering develops. The im¬ 
portance of magnetic resonances as a new spin relaxation 
mechanism, as pioneered in Ref. Ill] has been consolidated 
and further quantihed using extended models of more delo¬ 
calized magnetism, described within the Hubbard Hamilto¬ 
nian in the mean held approximation. Using efficient real 
space order N wavepacket propagation methods, spin re¬ 
laxation times in the nanosecond regime were obtained for 



FIG. 8: Main frame: Energy-dependence of spin (solid lines) and 
rescaled momentum (dashed line) relaxation times obtained from 
the real space wavepacket propagation method. Inset: momentum 
relaxation time derived the from single impurity model {I|] (solid 
line) and from the diffusion coefficient behavior (dashed green 
line). 


1 ppm hydrogen impurities, and for momentum scattering 
time few orders of magnitude shorter. More work is needed 
to describe larger hydrogen densities for which magnetic in¬ 
teraction and long range magnetic ground states could de¬ 
velop and interfere with relaxation effects. 


ACKNOWLEDGEMENTS 


The research leading to these results has received funding 
from the European Union Seventh Eramework Programme 
under grant agreement number 604391 Graphene Elagship. 
This work was also funded by Spanish Ministry of Economy 
and Competitiveness under contracts MAT2012-33911, the 
Severo Ochoa Program (MINECO SEV-2013-0295). S. R. 
Thanks the Secretaria de Universidades e Investigacion del 
Departamento de Economfa y Conocimiento de la Generali- 
dad de Cataluna. J.C.C. acknowledges the National Eund for 
Scientihc Research [E.R.S.-ENRS] of Belgium for hnancial 
support and the Research Concerted Action on Graphene 
Nanoelectromechanics (No. 11/16-037). Computational re¬ 
sources have been provided by the supercomputing facilities 
of the ’’Consortium des Equipements de Calcul Intensif ’ en 
Eederation Wallonie-Bruxelles (CECI), the UCL-CISM, and 
PRACE. J.E acknowledges support from DEG SEB 689 and 
GRK 1570. 


[1] D. Kochan, M. Gmitra, and J. Fabian, Phys. Rev. Lett. 112 , 
116602, (2014). 




















10 


[2] L. E. F. Foa Torres, S. Roche, and J. C. Charlier in In¬ 
troduction to Graphene-Based Nanomaterials: From Elec¬ 
tronic Structure to Quantum Transport (Cambridge Univer¬ 
sity Press, Cambridge, 2014). 

[3] K. P. Lob, Q.L. Bao, P.K. Ang, and J.X. Yang, J. Mater. Chem. 
20, 2277 (2010). 

[4] A.V. Krasheninnikov and F. Banhart, Nature Mater. 6, 723 
(2007). 

[5] J. Yang, M. Ma, L. Li, Y. Zhang, W. Huang and X. Dong, 
Nanoscale 6, 13301-13313 (2014). 

[6] D. C. Elias, R. R. Nair, T. M. G. Mohiuddin, S. V. Moro¬ 
zov, P. Blake, M.P. Halsall, A.C. Ferrari, D.W. Boukhvalov, 
M.l. Katsnelson, A.K. Geim, K.S. Novoselov, Science 323 , 
610 (2009); A. Bostwick, J.L. McChesney, K.V. Emtsev, T. 
Seyller, K. Horn, S.D. Kevan, and E. Rotenberg, Phys. Rev. 
Lett. 103 , 056404 (2009); J. Katoch et al., Phys. Rev. B 82 , 
081417(R) (2010); M. Jaiswal, C. Haley YX. Lim, Q. Bao, 
Ch.T. Toh, K. P. Lob, B. Ozyilmaz, ACS Nano 5 , 888 (2011); 
R. Jayasingha, A. Sherehiy, S.-Yu Wu, and G. U. Sumanasek- 
era. Nano Lett. 13 , 5098 (2013). 

[7] S. H. Cheng, K. Zou, F. Okino, H.R. Gutienez, A. Gupta, N. 
Shen, PC. Eklund, J.O. Sofo, J. Zhu, Phys. Rev. B 81 , 205435 
(2010); F. Withers, M. Dubois, and A. K. Savchenko, Phys. 
Rev. B 82 , 073403 (2010). 

[8] J. Moser, H. Tao, S. Roche, F. Alsina, C.M. Sotomayor Torres, 
A. Bachtold, Phys. Rev. B 81 , 205445 (2010). 

[9] S. Roche, D. Mayou, Phys. Rev. Lett. 79 , 2518 (1997); S. 
Roche, Physical Review B 59 , 2284 (1999); S. Roche, N. 
Leconte, F. Ortmann, A. Lherbier, D. Soriano, J.-C. Charlier, 
Solid. Stat. Comm. 152 , 1404 (2012). 

[10] J. Guillemette, S. S. Sabri, Binxin Wu, K. Bennaceur, P. E. 
Gaskell, M. Savard, P. L. Levesque, F. Mahvash, A. Guer- 
moune, M. Siaj, R. Martel, T. Szkopek, and G. Gervais, Phys. 
Rev. Lett. 110 , 176801 (2013). 

[11] O.V. Yazyev and L. Helm, Phys. Rev. B 75 , 125408 (2007); 
O.V. Yazyev, Phys. Rev. Lett. 101 , 037203 (2008); ibidem. 
Rep. Prog. Phys. 73 , 056501 (2010). 

[12] D. Soriano, N. Leconte, P. Ordejon, J.-Ch. Charlier, J.-J. Pala¬ 
cios, and S. Roche, Phys. Rev. Lett. 107 , 016602 (2011); N. 
Leconte, D. Soriano, S. Roche, P. Ordejon, J.-Ch. Charlier and 
J. J. Palacios ACS Nano 5, 3987 (2011). 

[13] S. Hemmatiyan, M. Polini, A. Abanov, A. H. MacDonald, and 
J. Sinova, Phys. Rev. B 90, 035433 (2014). 

[14] N. Tombros, C. Jozsa, M. Popinciuc, H.T Jonkman, B.J. Van 
Wees, Nature (London), 448 , 571 (2007). 

[15] W. Han and R. K. Kawakami Phys. Rev. Lett. 107 , 047207 
( 2011 ). 

[16] M. H. D. Guimaraes, A. Veligura, P. J. Zomer, T. Maassen, 1. 
J. Vera-Marun, N. Tombros, B. J van Wees, Nano Letters 12 
(7), 3512-3517 (2012) 

[17] T. Maassen, J.J. van den Berg, N. IJbema, F. Fromm, T. 
Seyller, R. Yakimova, and B. J. van Wees, Nano Letters 12, 
1498-1502 (2012). 

[18] B. Dlubak, M.-B. Martin, C. Deranlot, B. Server, S. Xavier, R. 
Mattana, M. Sprinkle, C. Berger, W. A De Heer, F. Petroff, A. 


Anane, P. Seneor, A. Fert, Nature Physics 8, 557-561 (2012). 

[19] P. J. Zomer, M. H. D. Guimaraes, N. Tombros, and B. J. van 
Wees, Phys. Rev. B 86, 161416(R) (2012). M. Drdgeler, F. 
Volmer, M. Wolter, B. Terres, K. Watanabe, T. Taniguchi, G. 
Guntherodt, C. Stampfer and B. Beschoten, Nano Lett 14 , 
6050 (2014). 

[20] 1. Neumann, J. Van de Vondel, G. Bridoux, M. V. Costache, F. 
Alzina, C. M. Sotomayor Torres and S. O. Valenzuela, Small 
9, 156160 (2013). 

[21] M.H.D. Guimaraes, P.J. Zomer, J. Ingla-Aynes, J.C. Brant, N. 
Tombros, and B.J. van Wees, Phys. Rev. Lett. 113 , 086602 
(2014) 

[22] S. Roche and S.O. Valenzuela, Journal of Physics D; Applied 
Physics 47 , 094011 (2014). 

[23] W. Han, R. K. Kawakami, M. Gmitra, J. Fabian, Nature Nan¬ 
otechnology 9, 794 (2014). 

[24] J. Fabian, A. Matos-Abiaguea, C. Ertler, P. Stano, 1. Zutic, 
Acta Physica Slovaca 57 565 907 (2007) 

[25] K. M. McCreary, A. G. Swartz, W. Han, J. Fabian, and R. K. 
Kawakami, Phys. Rev. Lett. 109, 186604 (2012). 

[26] M. Wojtaszek, 1. J. Vera-Marun, T. Maassen, B. J. van Wees, 
Phys. Rev. B 87 , 081402(R) (2013). 

[27] M. 1. Dyakonov and V. 1. Perel, Sov. Phys. Solid State 13, 
3023 (1971). 

[28] J. Balakrishnan, G. K. W. Koon, M. Jaiswal, A. H. Castro 
Neto, B. Ozyilmaz, Nature Physics 9, 284 (2013). 

[29] M. Gmitra, D. Kochan, J. Fabian, Phys. Rev. Lett. 110 , 
246602 (2013). 

[30] D. Van Tuan, F. Ortmann, D. Soriano, S. O. Valenzuela and S. 
Roche, Nature Physics 10, 857863 (2014). 

[31] G. Kresse and J. Furthmuller, Comput. Mat. Sci., 6, 15 (1996). 

[32] P. E. Blochl, Phys. Rev. B, 50 17953, (1994). 

[33] J. P. Perdew, K. Burke, and M. Emzerhof, Phys. Rev. Lett., 77 
3865, (1996). 

[34] R. J. Elliott, Phys. Rev. 96, 266 (1954).; Y. Yafet, in Solid 
State Physics 14 , ed. F. Seitz and Turnbull, p. 2 (1963); J. 
Fabian and S. DasSarma, Phys. Rev. Lett. 81 , 5624 (1998). 

[35] M. M. Ugeda, 1. Brihuega, F. Guinea, J. M. Gomez- 
Rodn'guez, Phys. Rev. Lett. 104 , 096804 (2010). 

[36] N.M. Peres, F. Guinea, A.H. Castro Neto, Phys. Rev. B 73 , 
125411 (2006); VM. Pereira, J.M.B. Lopes dos Santos, A.H. 
Castro Neto, Phys. Rev. B 77 , 115109 (2008). 

[37] Y.-W. Son, M.L. Cohen, S.G. Louie, Nature 444 , 347, (2006); 
J.J. Palacios, J. Femandez-Rossier, L. Brey, Phys. Rev. B 77, 
195428 (2008). 

[38] D. Soriano, F. Munoz-Rojas, J. Fernandez-Rossier, J.J. Pala¬ 
cios, Phys. Rev. B 81, 165409 (2010). 

[39] E.H. Lieb, Phys. Rev. Lett. 62, 1201 (1989). 

[40] M. Inui, S. A. Trugman, E. Abrahams, Phys. Rev. B 49, 3190 
(1994). 

[41] H. Ochoa, A. H. Castro Neto, F. Guinea, Phys. Rev. Lett. 108 , 
206808 (2012). 

[42] D. Kochan, M. Gmitra, and J. Fabian, arXiv: 1306.0230 


